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A time-dependent variational Lagrangian analysis based on the Gross-Pitaevskii energy func- 
■ tional serves to study the dynamics of a metastable giant vortex in a rapidly rotating Bose-Einstein 

' condensate. The resulting oscillation frequencies of the core radius reproduce the trends seen in 

' recent experiments [Engels et ai, Phys. Rev. Lett. 90, 170405 (2003)], but the theoretical values 

, are smaller by a factor ~ 0.6-0.8. 
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^ . One of the most characteristic features of neutral superfluids is their inability to rotate uniformly with solid-body 
velocity Vgh = Jl x r. Instead, a rotating superfluid contains one or more singly quantized vortex lines with circulation 
, h/M , where M is the atomic mass. This phenomenon is familiar in rotating superfluid ^He 1], although visualization 
of the resulting vortex arrays is quite difficult The creation of the first Bose-Einstein condensate (BEG) in dilute 
^^Rb gas 13 challenged experimentalists to demonstrate the existence of similar quantized vortices in dilute BECs 
■ O ' and to study their properties. After the development of techniques for creating a single vortex f3| and small vortex 
C ' arrays j5| , rapid progress led to large vortex lattices and rotating condensates with large angular momentum jjj ■ 
One recent method to increase the angular momentum has used an intense focused laser beam to remove the 
central region of a rotating condensate 'sl^^'o'l. This method notes that the particles near the center have relatively 
little angular momentum, so that the resulting annular condensate has an increased angular momentum per particle. 
In this way, it has been feasible to reach rotation rates f2 « 0.97 wj,, where wj, is the radial confining trap frequency. 
This process creates a relatively long-lived "giant vortex core" that contains a large number of phase singularities, 
resulting in a macroscopic circulation around its boundary. The remaining annular fluid has a dense vortex lattice 
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I. INTRODUCTION 
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that rotates at the angular velocity fl. In addition, the radius of the core oscillates at a frequency of order 3-3.5 uj_i_, 
^\ , with relatively weak damping |9| . Recently, a theoretical study of this behavior has relied on numerical simulation 
' of the time-dependent Gross-Pitaevskii (GP) equation in two dimensions ^3 ^^'^ found an oscillation frequency 
~ 2.6 w_L, somewhat smaller than the observed values. The present work provides an alternative theoretical study 
' of the same phenomenon using a variational Lagrangian that contains the inner radius and the amplitude of the 
radial velocity as time-dependent parameters. The resulting Lagrange's equations yield an undarnped oscillation with 
definite amplitude and a frequency that varies between 0.6 and 0.8 of various measured values [MS- Although the 
' two theoretical approaches appear quite different, they both rely on the GP equation and both predict oscillation 
■ frequencies that are definitely smaller than the observed ones. 
^ ' Section UTI reviews the variational Lagrangian procedure and applies it to the present problem. Sections IIIII and Hvl 
O , contain analyses of two specific models. The first is a quasi-two-dimensional Thomas-Fermi (TF) condensate with 
O tight confinement in the z direction; the two-dimensional character of the condensate wave function simplifies the 
analysis considerably. The second is a more realistic three-dimensional TF condensate, but the final results are rather 
similar in both cases. The comparison with the experimental observations is discussed in Sec. IVl 



II. BASIC PROCEDURE 



The approach here is the variational Lagrangian method, in which a trial condensate wave function \1/ includes 
various parameters whose time dependence follows from Lagrange's equations for the Lagrangian 

L^-m dV U*—^^.— ]-Em, (1) 
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where E[^] is the energyfunctional. This flexible procedure has served successfully to study the low- lying collective 
modes of a condensate (Tj.] and the precession of a single vortex line 0, 0, 0, ^( in a trapped condensate. It 
is convenient here to work in the laboratory frame of reference, and to impose conservation of particle number and 
angular momentum explicitly. 

The experiment 8] studies N ^ 10^ atoms of ^^Rb in an axisymmetric trap with uj±/2tt = 8.3 Hz and u!z/2tt — 5.4 
Hz, so that the relevant oscillator lengths are d± ~ 3.73 /im and dz ~ 4.62 /im. Thus the nonrotating condensate 
is somewhat elongated, with TF dimensions Rz{0) / R±{0) = uj±/uJz ~ 1.54. When the condensate rotates rapidly, 
however, it expands radially and contracts axially according to 16] 

wheieU^n/uj^. Taking = 0.97 0, the standard TF expressions ^ yield R±m/d± w 19.9 and w 7.45, 

with Rz(Ti)/Rj_(n) = {l-Tfy/^Rz{0)/R±{0) w 0.37. Thus the rotating condensate is definitely disk shaped. 
In the TF limit, the gradient of the density j^Pp is small, and the energy functional becomes 

E[^]^ JdV (^lMv'\^\' + Vu.m'+^-^\m'y (3) 

where v'^ is the squared velocity of the superfluid, Vtr is the trap potential, and a « 5.77 nm is the s-wave scattering 
length. The density j^'p is taken as the familiar quadratic TF distribution with an empty core of radius Tc 

1*1' = |*TFp^(r-re). (4) 

The core radius Tc serves as one dynamical variable, and the resulting radial motion requires a time-dependent radial 
velocity Vr = {h/M) dS/dr in addition to the azimuthal component — (h/AIr) dS/d(t), where S is the phase of the 
condensate wave function. 

III. QUASI-TWO-DIMENSIONAL CONDENSATE 

To simplify the analysis, the condensate is first assumed to be thin in the axial direction with a Gaussian dependence, 
so that '^{r,z) « 5'2d(i") V'ol^)- Here, ^o{z) = {y^dz)~^^^ exp{—^z'^ /dl), where dz = ^JhlMu^z- For the two- 
dimensional wave function \I'2d(r), it is convenient to use dimensionless variables with frequencies (and inverse time) 
scaled by the radial trap frequency w^, energies scaled by fiwj^, distances scaled by the radial oscillator length 
d^ — \Jh/Mijj^_, and the condensate wave function ip = ^'/-s/iV normalized to 1. 

The trial wave function used here corresponds to a radial TF two-dimensional density for r < R with the central 
region of radius Tc removed, so that 

IV'zdMP = (^1 - ^) 0iR - r) 0{r - r,), (5) 

where R = R±/d± is now dimensionless and C2d is a normalization constant. The phase S of the wave function is 
assumed to have the radial part 

Sr^a{^r^~rR) (6) 

to allow for the oscillatory radial motion, where a is a time-dependent parameter. In addition, the phase has the 
familiar contribution Sy from the vortices plus the circulation around the hole representing the missing vortices 

S'„(r) =^arctan(^^— ^) +i/o(/), (7) 

j V ^ ^3 / 

where the sum is over all vortices j. 

The radial velocity Vr = dSr/dr = a (r — R) vanishes at the outer boundary, which remains fixed in this model. 
For simplicity, only Tc and a are taken as dynamical variables. In contrast, the vortices are not treated as dynamical 
variables; instead, they are assumed to have uniform dimensionless areal density p and to move with the local azimuthal 
velocity. In detail, the azimuthal component of the dimensionless velocity is 
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where vq = Trpor^ is the constant quantum number of the circulation around the hole of initial radius tq with initial 
vortex density pg. The sum can be approximated by an integral 



P 



(fr'- 



r — r cos ( 



j,2 _ 2rr' cos 0' + r'- 

2 



r 



2ttp 
r 



r 



'^c , I'd 

npr — np \ 

r r 



(9) 



Note that the initial azimuthal component of velocity irpor represents a pure solid-body rotation at a dimensionless 
angular velocity Q, = irpo since the last two terms of Eq. cancel for Tc = rg (when p = po). 

The analysis is carried out in the laboratory frame, and the conservation of particle number and angular momentum 
are explicitly incorporated into the resulting Lagrangian. It is straightforward to obtain the normalization constant 



7ri?2 (1 _ 2/2)2 ' 



(10) 



where y = rd R determines the instantaneous position of the core radius. Similarly, the angular momentum per 
particle (in units of H) is 



hd^ I Srrv^\'<ij2d\^ ^\npR'{l-y^) + vo = \T^PQR'{'^ + 'iyl). 



(11) 



expressed in terms of the initial parameters po and y^. Since and are taken as constants, this relation determines 
the actual vortex density 



3 hd — ^0 



ttR^ l-y2 



(12) 



as the core radius y changes with time. Note that the total number of vortices 7rpR^{l — y"^) — TTp{R? — r^) in the 
annular region Vc < r < R remains constant, so that the vortex density changes only because the inner radius Tc 
changes. Correspondingly, the initial conditions fix the quantity l2d — = ^R^Vl{\ — y^). Assuming typical values 
i? « 20, n « 0.97, and yo « 1/3, I find hd ~ 158 and i^o « 43. 

After the z integration is carried out with the tightly confined Gaussian, the dimensionless Lagrangian per particle 
becomes 



L2d = \i 



'^2d^, V'2d 



dt 



dt 



E2d[ip: 
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where 



E2dW 



2d 



\{vl + vl) |^2dP + \ A^2d? + ^/2^^ |^2d| 



(14) 



is the two-dimensional energy functional. All the integrals can be performed analytically and Eqs. (|10|) and (|ll|l ensure 
the conservation of particles and angular momentum. Note that the same Lagrangian also describes a condensate 
that is uniform along a length L in the z direction. In this case, the coefficient of the quartic term in i?2d is changed 
to 2'KNa/ where N /L is the number of particles per unit length. 

The resulting Lagrangian depends on two dynamical variables: a that fixes the radial component of velocity and 
y = rd R that fixes the position of the inner core radius. It is convenient to use the equivalent quantity 



C2d = -i2d/i?' - «[/2d(y) - 2] + f2d{y) + £2d{y) 



(15) 



where the first term arises from the explicit time derivative in 1)13(1 . the second is the radial kinetic energy, and E2d 
constitutes the remaining terms in the energy. It has the explicit form 



£2d(y) = n {l-ylf 



52d(y) 
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h{yf32d{v) 



1. „ 2^ 4V2 Na 1 
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where the first term (proportional to ) is the azimuthal kinetic energy, the second is the trap energy, and the last 
is the interaction energy (note that this term is typically small, despite the large value of Na/dz ~ a few thousand, 
owing to the factor ~ 6 x 10^^). To be very specific, Cid involves the following functions 

^''^y^ 15(1 + ' (^^^ 



IV. THREE-DIMENSIONAL THOMAS-FERMI CONDENSATE 



Before studying Lagrange's equations for C2d, it is convenient to consider the corresponding case of a three- 
dimensional Thomas-Fermi condensate, where 

»rHr.^)f . CJ, (l - ^ - Ij) . (l - ^ - Ij) - .J. ,21, 

using the same scale factors lu±_ and dj_ to define dimcnsionlcss variables. The normalization integral gives 

" 8^i?2^,(l-y2)5/2' (22) 
which can be compared to Eq. (|10|l . Similarly, the angular momentum becomes 

Itf ^^7TpR\l-y^) + ,^o, (23) 
and the conservation of angular momentum then yields the vortex density 

7 /tf^ (24) 

All the terms in the Lagrangian can be evaluated analytically, and the resulting expression has the same form as 
Eq. ifT^ but with different functions 

Ctf = a[fTF{y) - ^] + a^frpiy) + Sxriy)- (25) 



In particular, 



and 



^ 9 + 5y2 5 7r+|2/(l-y2)i/2(3_i4y2 + 8y4)_2arcsiny 

/^Hy) = ^^ ^ O^TjHT^ (26) 



fTF(2/) = 1^ (1-yo')' 



9TF{y) + -f-^ + KvYiTFiy) 
1 - y 



-[2 + 5,^ + (l-f7 )(l-,^)] + --^^^-^, (27) 

The extra factor of in the interaction term will be seen to eliminate the dependence of the number of atoms 
(assuming only that the condensate remains in the TF limit). In addition, the remaining functions have the altered 
forms 

2 + 5y^ 

9TF{y) = ^^(^_^2)2 ^ (28) 

5 (1 ~ y^f'^ (-4 + y^) - 31ny -f 31n[l + (1 - y^f'^] 
^^^^y^ ^ 6 (1-^2)5/2 • (29) 
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V. DISCUSSION OF RESULTING MOTION 



(33) 



Variation of the effective Lagrangian L with respect to a yields the dynamical equation 

2) ^-2/(2/) a, (30) 
dy 

which expresses a in terms of y and y. Similarly, variation with respect to y yields a more complicated equation 

a^ = (" + «)^ + ^-0 (31) 

that expresses a and a in terms of y. 

Unfortunately, it is not easy to interpret this pair of coupled dynamical equations, and it is preferable to assume 
small oscillations about the initial conditions j/(0) = yo a-nd j)(0) = (namely the hole in the rapidly rotating 
condensate starts from rest). Consequently, a itself is a small quantity, as is 5 — y^ — y, which is positive when the 
hole contracts from its initial scaled radius yo- In this way, Eq. H3U|) becomes 

' = "7(^"' ^''^ 

where the prime denotes a derivative with respect to y. Similarly, Eq. (|31|l has the linearized form 

. ^ -E'{y^)^5£"{y^) 

nyo) 

A combination of these two linearized equations yields a simple dynamical equation 

6 + uj^6 = Lj^A, (34) 

where 

2 2£"{yo)f{yo) 
uj = 5 — , (35) 

Wiyo)]' 

and 

They depend on the initial radius yo = rg/R and on the small parameters Na/dzR^ in the two-dimensional case and 
Na/ d±RzR^ for the three-dimensional TF case (these latter contributions come from the interaction energy, which is 
small because the expanded outer radius reduces the particle density far below its value for a nonrotating condensate) . 
Equation (|34|l has the solution 

5{t) ^ A{1- cos ut). (37) 

It represents an inward oscillation of the core radius with angular frequency lj (measured in units of lo^_) and amplitude 
A (like y, A is measured in units of i? = R±/d±). As a result, conservation of angular momentum and the parabolic 
trap potential dominate the dynamical motion because the interaction term makes only a small contribution; this is 
the effect of the Coriolis force mentioned in In all cases, the amplitude A is relatively small, partially justifying 
the expansion through quadratic contributions to C. Note, however, that the minimum core size is tq — 2AR, which 
corresponds to a significant shrinkage for the smaller rotation rates Q < 0.9. 

Table I contains typical calculated values for some of the experimental parameters used in p| , using the functions for 
a quasi-two-dimensional condensate. Since the present theory is rather crude, however, the general three-dimensional 
TF radii have been used even for the rotating quasi-two-dimensional condensate. The choice N — 2.2 x 10^ corresponds 
to Fig. 4 of Ref. 0, which shows a rotating condensate with = 0.9; the core executes weakly damped oscillations 
with dimensionless frequency uj/u!± = 3.5. This frequency is higher than the value co/cu^ = 2.32 found here for a core 
radius 1/3 of the condensate radius, but the dependence of the oscillation frequency and amplitude on the rotation 
frequency f2 both show the observed trend (namely, "decreasing the initial rotation rate of the condensate leads to 
faster core-oscillation frequencies and increases the amplitude of this oscillation" (sj). The corresponding values for 
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N = 3.0 X 10^ illustrate the extremely weak dependence on the number of atoms. In fact, if strict two-dimensional 
TF radii were used, the dependence on N would cancel completely, as will be seen to hold for the three-dimensional 
TF condensate. 

Table II is similar to Table I except that it is evaluated for the three-dimensional Thomas-Fermi condensate density. 
In this case, it is not difficult to use the standard TF expressions for the condensate radii JJJ to verify that the ratio 
{15Na/d±){R'^Rz)^^ in fact reduces to the value 1 — fi, so that the N dependence indeed drops out. In all the cases 
considered, the three-dimensional TF condensate yields somewhat lower oscillation frequencies and smaller amplitude 
for the motion. 

To provide a more detailed comparison between the theoretical predictions and experimental values. Table III 
contains 5 sets of measured frequencies along with the corresponding predictions of both models. For the quasi- 
two-dimensional condensate, the ratio of predicted value wth to measured value Wobs varies between 0.61 and 0.79, 
whereas for the three-dimensional TF condensate, the same ratio varies between 0.49 and 0.76. In both cases, the 
agreement improves considerably at the fastest rotations. This trend may reflect the present use of a spatially averaged 
velocity for that induced by the vortices and the macroscopic circulation, ignoring the discrete character of the vortices; 
such an approximation should work best at large f2. 

One possible explanation for the discrepancy between the experimental values and the results of the present vari- 
ational Lagrangian formalism is that the actual dynamical condensate does not conserve the number of vortices; in 
this case, angular momentum would not strictly be conserved. At present, it is not clear how to incorporate such 
additional dynamical variables into the analysis. In any case, it is notable that the numerical simulation of Simula 
et al. ^3 also seems to predict an oscillation frequency that is too small by a comparable amount. Conceivably, the 
quite dramatic perturbation caused by the intense laser beam in burning a relatively large hole may well produce 
significant depletion, so that a theory based solely on the time-dependent GP equation would need to be modified to 
include noncondensate contributions. 
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TABLE I: Calculated parameters for a quasi-two-dimensional condensate 



iV = 2.2 X 10^ iV = 3.0 X 10'' 



r-o/R 




R = R±/d± 


Lj/uj± 


A 




A 


1/3 


0.90 


13.20 


2.32 


0.097 


2.34 


0.099 


1/3 


0.94 


15.29 


2.19 


0.063 


2.20 


0.064 


1/3 


0.97 


18.74 


2.10 


0.032 


2.10 


0.033 


1/5 


0.90 


13.20 


2.50 


0.081 


2.52 


0.083 


1/5 


0.94 


15.29 


2.32 


0.054 


2.33 


0.056 


1/5 


0.97 


18.74 


2.18 


0.029 


2.19 


0.030 



TABLE II: Calculated parameters for a three-dimensional TF condensate (the predicted frequencies and amplitudes are inde- 
pendent of A'' as long as the TF approximation remains valid) 



ro/R = 1/3 



ro/R = 1/5 



Q = Sl/tiJx 

0.90 
0.94 
0.97 



R = R±/d^ 

13.20 
15.29 
18.74 



A 




2.20 
2.15 
2.11 



A 

0.052 
0.034 
0.018 



TABLE III: Comparison of selected experimental values with theoretical predictions for two-dimensional and three-dimensional 
models 



experimental values" 



quasi-two-dimensional model 



three-dimensional TF model 



ro/R'' n = 


= Q/u)± 


N 


OJobs/l^x'^ 






A 




I^th/<^obs 


A 


0.40 


0.79 


2.8 X 10*^ 


4.5 


2.74 


0.61 


0.171 


2.17 


0.49 


0.122 


0.24 


0.88 


2.3 X 10^ 


4.0 


2.48 


0.62 


0.098 


2.16 


0.54 


0.064 


0.32 


0.90 


2.1 X 10*^ 


3.5 


2.32 


0.67 


0.095 


2.09 


0.60 


0.060 


0.32 


0.95 


2.5 X 10*^ 


2.9 


2.17 


0.75 


0.053 


2.06 


0.71 


0.032 


0.35 


0.96 


2.3 X 10® 


2.7 


2.13 


0.79 


0.044 


2.05 


0.76 


0.027 


"Reference 9 
''Reference j£ 
■^Reference 


assigns 
1 assigns 


an uncertainty of ^ 20% to the initial core radius ro. 
an uncertainty of 6 Hz to the first value of uiahs (37 Hz), namely ~ 


16%. In contrast, the other 


experimental 



oscillation frequencies have uncertainties of < 1 Hz, which is < 4% even for the last value. 



